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Abstr  acrfc 

Using  the  basic  principles  of  operation  of  the  ter- 
restrial hydromagnetlc  dynamo.,  it  is  shown  that  the  solar  convec- 
tive zone  will  generate  traveling  dynamo  waves  consisting  of 
toroidal  and  poloidal  magnetic  fields*  The  waves  migrate  from 
the  poles  toward  the  equator  near  the  top  cf  the  connective  zone 
and  in  the  opposite  direction  near  the  bottom-  The  thinness 
of  the  convective  zone  is  responsible  for  a migratory  dynamo 
in  the  sun.f*ather  than  a stationary  dynamo  as  in  the  Earth.  An 
investigation  of  the  dynamics  of  hydromagnetlc  dynamos  indicates 
that  the  poloidal  field  will  be  strongest  above  middle  latitudes 
and  the  toroidal  field  below.,  in  agreement  with  the  observed 
nonuniform  rotation  and  the  appearance  of  sunspots  in  low 
latitudes. 


1.  Introduction 


The  migration  toward  the  equator  and  the  reversal  of 
polarity  of  sunspots  each  eleven  year-  half -cycle  are  probably 
the  most  conspicuous  aspects  of  the  overall  solar  magnetic  acti- 
vity. The  emission  of  excessive  radio  noise  and  charged  particles, 
prominence  and  flare  activity,  and  coronal  streamers,  as  well  as 
direct  observation  of  the  magnetic  field  over  the  disk  of  the  sra, 
(Babcock,  1953)  indicate  that  the  solar  magnetic  field  is  complex  - 
certainly  not  just  a dipole  - and  has  a nearly  periodic  time 
dependence. 

The  Earth's  field, which  is  predominantly  a dipole,  with 
its  moderate  fluctuations  in  space  and  time,  is  satisfactorily 
accounted  for  by  a stationary  hydrcmagne tic  dynamo  model  based  on 
the  convective  motions  in  4'he  conducting  fluid  core  and  the 
Coriolis  forces  resulting  from  the  Earth's  diurnal  rotation 
( Elsasser,  1946/7 } Parker,  1954).  The  two  necessary  dynamo  con- 
ditions, viz.  rotation  of,  and  convection  within,  a large  elect- 
trically  conducting  fluid  body,  are  satisfied  by  the  sun.  This 
suggests  that  the  solar  magnetic  field  is  the  result  of  a hydro- 
magnetic  dynamo*  it  will  be  our  purpose  here  to  construct  a 
crude,  appropriate  model. 

We  have  at.  our  disposal  a convective  layer  near  the  surface 
of  the  sun  with  a thickness  of  the  order  of  10  km.  The  convection 
results  in  the  first  place  in  nonuniform  rotation  with  the  outer 
parts  of  the  convective  layer  lagging  behind  the  inner.  These 
fluid  motions  also  seem  to  produce  a poloidal  magnetic  field 
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which  ha-  only  a small  dipole  component,  which  is  of  opposite 
sign  in  the  northern  and  southern  hemispheres,  which  undergoes 
continual  migration  toward  the  equator-  from  both  poles,  and  which 
reverses  sign  at  any  given  point  on  tbe  sun  every  eleven  years. 
This  suggests  ’’waves”  of  magnetic  flux  originating  near  the  poles 
and  being  propagated  regeneratively  toward  the  equator  where  they 
are  dissipated;  the  period  of  the  waves  would  be  22  years. 


2.  Qualitative  Discussion 


Before  developing  an  elaborate  analysis  of  a dynamo  model 
satisfying  these  requirements , let  us  develop  the  model  in  a 
qualitative  way  in  order  to  understand  physically  the  hydromagnetic 
processes  we  might  expect  to  occur.  (The  calculations  of  the  next 
section  will  be  of  interest  to  the  person  interested  especially 
in  dynamos.)  Figure  1 shows  a slab  of  fluid  representing  a 
piece  of  the  convective  shell  of  the  sun.  At  time  t we  take 
the  sectioq  to  be  a rectangular  parallelepiped,  neglecting  the 
curvature,  as  indicated  by  the  broken  lines.  At  some  later  time 
the  nonuniform  rotation  will  have  sheared  the  section  into  the 
form  shown  by  the  solid  lines,  and  a loop  of  flux  initially  in 
a meridional  plane  will  be  carried  by  the  fluid  into  the  nearly 
horizontal  position  represented  by  the  unbroken  elliptical  band. 
This  is  just  the  familiar  process  of  the  nonuniform  rotation 
generating  a toroidal  field  from  the  poloidal  field.  The  net 
pclcidal  field  is  unchanged  by  the  process. 

Consider  next  a meridional  section  through  the  convective 
zone  illustrated  in  figure  2(a).  Two  initial  magnetic  loops  in 
the  meridional  plane  are  represented  by  the  heavy  closed  lines. 

We  regard  the  meridional  section  illustrated  in  figure  2 as  the 
left  hand  end  of  the  slab  in  figure  1,  so  that  the  nonuniform 
rotation  displaces  the  top  side  of  the  meridional  section  in 

figure  2 out  of  the  paper  relative  to  the  bottom.  The  toroidal 
field  will  be  out  of  the  paper  near  the  center  and  into  the  paper 
at  the  right  and  3oft  hand  ends,  as  indicated  by  the  plus  and  minus 
signs  respectively.  The  convective  motions  interacting  with  the 
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toroidal  fields  will  produce  magnetic  loops  with  nonvanishing 
projections  on  the  meridional  planes  (Parker,  1954).  The  sense 
of  these  loops  depends  upon  the  sense  of  the  toroidal  field.  Let 
us  assume  that  when  the  toroidal  f ield  is  directed  out  of  the  paper 
the  individual  loops  are  in  a counterclockwise  direction)  when 
the  toroidal  field  is  into  the  paper,  the  loops  will  be  in  a clock- 
wise direction.  Many  such  loops  will  be  produced,  and  it  can  be 
shown  that  they  will  coalesce  into  large  loops,  which  are  indi- 
cated by  the  light  lines  in  figure  2(a).  If  now  we  add  the  newly 
formed  large  loops  to  the  initial  loops,  we  see  from  the  sense  of 
rotation  that  the  initial  loops  will  be  reinforced  in  the  vicinity 
of  B and  D and  degenerated  at  A and  C.  Figure  2(b)  shows  the 
result  of  this  superposition,  resulting  in  the  initial  loops  being 
displaced  to  the  right.  Thus,  we  have  a dynamo  which  not  only  may 
have  regenerative  feedback  between  the  poloidal  and  toroidal  fields, 
but  which  also  produces  a migration  of  the  field*  The  direction 
of  the  migration  depends,  of  course,  on  the  sense  of  the  loops 
and  of  the  toroidal  field. 


3,  The  Migratory  Dynamo  Equations 


Consider  the  space  determined  by  the  rectangular  coordin- 
ates ( x,  y,  z ) and  filled  with  a fluid  with  magnetic  viscosity 

/ 

We  shall  assume  that  the  fluid  has  several  motions:  There  is  a 

motion  eKv(x.,H)  , where  , e,  , and  g,  are  unit  vectors  in  the 
x,  y and  z directions  respectively.  3 will  represent 
the  shearing  of  the  nonuniform  rotation  found  in  stellar  convec- 
tive zones*  Besides  v(x,z)  we  will  assume  that  the  fluid  is 
convecting  in  columns  parallel  to  the  ? -axis.  The  columns  are 
rotatixig  about  thsir  axes  as  they  would  in  an  actual  star  because 
of  Coriolis  forces.  To  compare  this  idealized  rectangular  model 

l 

with  the  actual  physical  situation  in  a star,  we  will  regard  some 
cell,  say  O < =*■  < a j O < z as  a piece  of  the  convective 

shell  in  which  the  curvature  has  been  neglectedj  we  will  think 
of  the  positive  x-axis  as  pointing  to  the  south,  the  y-axis  to 
the  east,  and  the  z-axis  upward. 

We  shall  represent  the  toroidal  field  by  JB  , where 

fe  = Sr  (r,  t)  , r = ^ (1) 

p«l^idls3L  field  is  more  easily  represented  by  the  vector 
potential 

f\  *=  A (jTj  t)  (2) 

than  by  the  magnetic  field  itself,  which  has  both  z and  x com- 
ponents. We  shall  take 

isA 


SB 


= O 


(3) 
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The  interacting  with  the  toroidal  field  B of  a rotating 
convective  column  moving  in  the  z direction  generates  a magnetic 
loop  with  a nonvanishing  projection  on  the  xz- plane.  The  details 
of  the  loop  generation  may  be  found  in  an  earlier  paper  (1954). 
Briefly,  what  occurs  is  that  a rotating  or  cyclonic  convective 
column  rises  up  through  the  toroidal  field  poking  the  field  ahead 
of  it  to  make  a bump.  The  rotation  of  the  convective  column,  then 
twists  the  bump  through  90°  giving  magnetic  loops  in  the  meridional 
plane.  The  loops  coalesce  to  form  the  poloidal  field.  Since  we 
represent  the  poloidal  field  by  its  vector  potential,  each  loop 
will  appear  as  a hump  in  the  function  A ( JC , t)  /.  If  we  assume 
that  the  rate  of  generation  of  at  any  giver  point  is 

proportional  to  the  toroidal  field  from  which  it  is  generated,  we 
obtain 

?> & - * VvA  = Hr)  b (4) 

2>t  K 1 

where  P (£)  is  some  function  characteristic  of  the  convective 
motions. 


The  toroidal  field  is  generated  through  the  interaction  of 
v(x,z)  with  the  poloidal  field.  Hence 

- t7*  [X  * (?X  A)]  (5) 

But  since  and  have  x and  z components  only 


\7*  [*  x $7  «Aj]  - ~ 


(6) 


Hence 


2>t 


(yv)x(vAj 


/ Tt  \ 

V ' ) 
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(4)  and  (7)  represent  the  complete  dynamo  equations. ; To 
effect  their  solution  we  define  k (k)  and  ^(tOa3  the 

Fourier  transforms  of  A (ciO,  (jr)^  and  F (jf  ) respectively, 

so  that  ..c  l- 


A ( £,  i ) =•  j dk  /dka  exp6k-*J  a Cfe,  0 

— qd 


_ 4 oo  -fees  ^ 

QbA)  ~ J dx.  j d a ejp(-  { jk’jr)  A (.f  j t ) ! 


\ 


/ o \ 
’« / 


etc.  where 


k *=  W,  +•  jga 


(9) 


We 

note 

r r 

400 

1 

4nx 

M- 

1b  e.>.| 

-CO 

and 

-fco  -M* 

J-  Lx/ia  «xp(-;t4R)xW= Jdkfdlf  fc(fc-fc9^*'a(k') 

/-C  ''*• 


Operating  on  (4)  and  (7)  with 

. +-<*>  -+  o° 


-L.  cl. 


»7tj 


d B e^p^-  ih‘  jr) 


we  obtain  the  two  equations 


J=L  a (fc,  t)  4-  i;1  a (i,^)  = J dk,  / d Ic  i.  Cjfc  J)  Cjb~  "ty  v 10) 

-J-tm 


3.b(fc,t)  **Ob(|=,t)  - Adtl'  f dk>'  k(k-k)*:l  'aCfc',0  (11) 

"&t  >-«  y - «, 


#7  T!he  integral  eguation  (1-3)  was  first  derived  by  the  author 
from  a direct  consideration  of  the  toroidal  and  poloidal  fields 
and  their  interaction  with  the  fluid  velocity.  After  reviewing 
this  rather  lengthy  derivation.  Professor  W»  M.  El3as3er  pointed 
out  that  the  interaction  could  be  condensed  to  the  two  differential 
equations  (4)  and  (7),  thereby  greatly  simplifying  the  development. 
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The  general  solution  of  (10)  is 


a fei)  « f jk/  L:  r{k)  ±>(k-k',  o (i2) 

S—  ft)  - oo 


Putting  (12)  into  (11)  we  obtain  finally 


■a.  tCk.O  = 

1st 
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4.  Traveling  Dynamo  Waves 


Let  us  consider  the  solution  of  the  dynamo  equation,  (13). 

To  establish  that  traveling  dynamo  waves  exist,  consider  the  case 
where  the  convective  motions  and  shearing  of  the  nonuniform  rotation 
are  independent  of  position.  Then 


w = h , r<»  = r 

where  H and  P are  constants.  Then 

h (fe)  = £*  H §>  (h)  , « r 8 Qk) 

where  S(fc)  is  a Dirac  delta  function*  (13)  reduces  to 

.t 


(14) 


(15) 


"at 


^ k*  b (fc,  ( ; + i Hr  k,  f JtWpf'k'C  t'-tjjbfejl')  (16) 

~ -S~m 

It  is  readily  shown  that  the  solution  of  (16)  may  be 


written 


- pCk) 


l 


J 


*■)  (jk)  ex 


p <- 


(IV) 


The  + has  the  same  sign  as  k,  if  H P )>o  and  the  opposite  sign 
if  HP<  O • (IV)  represents  a plane  wave  oropagating  with 
velocity  [|H  P / (2k,)]  < in  the  Ik  direction.  The  second  term 
represents  degeneration  of  the  existing  field  and  vanishes  for 
large  values  of  t . Thu^  for  large  values  of  t , (IV)  reduces  to 

b^U-i-CO  .kjt  * <18> 

This  represents  a plane  wave  traveling  in  the  negative  Si  .direction 
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if  HT  >C  and  in  the  positive  jk  direction  if  H P < O . The 
amplitude  of  the  mode  with  wave  number  jk  is  constant  if 

I h r/  = 

k 

the  amplitude  grows  exponential  with  time  if  IHFl  is  greater 
than  k and  decays  if  it  is  less. 

From  (13)  we  see  in  general  that  the  modes  Jk.  and  jfk  are 
coupled  through  the  quantities  jb  (k  - i=/)  ~ jk*^  . Thus  the 

modes  will  not  be  independent  except  in  the  special  case  of  (14). 
Coupling  between  modes  increases  the  difficulty  of  solving  (13) 
and  the  general  case  is  discussed  only  briefly  in  the  appendix, 
where  two  special  cases  are  worked  out. 
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5.  The  Kinematics  of  the  Solar  Dynamo 

Consider  the  cyclonic  motion  of  a rising  convective 
column  in  the  northern  hemisphere  of  the  solar  convective  zone. 

The  Coriolis  forces  resulting  from  the  influx  of  matter  at  the 
bottom  end  of  the  column  give  rise  to  a counterclockwise  rotation 
about  a vertical  axis  as  viewed  from  outside  the  convective  zone. 
Unless  dissipated  by  viscosity  the  rotation  of  the  fluid  will 
continue  until  it  is  brought  to  a halt  by  the  Coriolis  forces  of 
opposite  sense  due  to  the  efflux  at  the  top  of  the  column.  This 
was  the  situation  assumed  in  the  core  of  the  Earth,  illustrated 
in  figure  3(a).  The  convective  motions  in  the  sun,  however,  are 
much  more  rapid  than  in  the  core  of  the  Earth,  1 m/sec  as  compared 
to  0.1  mm/ssc,  and  of  much  larger  scale,  with  the  result  that  the 
solar  Reynolds  numbers  are  large.  Thus,  the  motions  in  the  solar 
conxrective  zone  constitute  well  developed  turbulence  in  contrast 
to  the  irregular  convective,  though  not  fully  turbulent,  fluid 
motions  in  the  core  of  the  Earth.  The  cyclonic  motions  in  the  sun 
have  to  contend  with  eddy  viscosity  besides  the  relatively  neg- 
ligible molecular  viscosity.  Now  an  eddy  in  a fully  developed 
field  of  turbulence  has  a lifetime  of  the  order  of  its  diameter 
divided  by  its  velocity,  i.e.  an  eddy  lasts  about  l/x  of  one 
’’revolution" . Thus,  we  expect  ihe  cyclone  resulting  from  the  in- 
flux at  the  base  of  a rising  convective  column  in  the  sun  to  be 
dissipated  sometime  before  it  has  made  one  complete  revolution. 

The'  efflux  at  the  upper  end  of  the  column  will  then  give  rise  to 
a new  rotational  motion,  this  time  an  anticyclone.  This  means 


that  in  the  northern  hemisphere  of  the  solar  convective  zone  we 
expect  a counterc  lockwise  rotation  to  be  associated  with  an  influx 
of  matter  to  a convective  column  and  a clockwise  rotation  with  an 
efflux,  as  viewed  from  outside  the  sun.  Vice  versa  in  the  southern 
hemisphere.  The  rotation  of  the  fluid  in  a rising  convective 
column  in  the  northern  hemisphere  is  illustrated  in  figure  3(b) 

To  state  this  matter  a little  differently,  we  note  that 
an  influx  is  associated  with  the  lower  half  of  rising  columns, 
and  an  efflux  with  the  upper  half.  Hence,  in  the  outer  half  of 
the  solar  convective  zone,  rising  and  sinking  columns  are  assoc- 
iated with  clockwise  and  countgrclo  okwise  rotations  respectively, 
and  vice  versa  in  the  inner  half.  Figure  4 illustrates  the  loops 
that  will  be  produced  in  the  inner  and  outer  halves  of  the  con- 
vective zone  assuming  a toroidal  field  running  from  west  to  east. 

We  note  that  in  the  outer  half,  the  loops  due  to  both  sinking  and 
rising  columns  are  in  a counte  rolockwise  direction  and  in  the 
inner  half  clockwise,  as  viewed  from  the  east. 

Consider  the  hydromagnetic  dynamo  that  will  result  from 
the  convective  motions  outlined  above.  As  mentioned  earlier  the 
turbulence  within  the  convective  zone  gives  a large  magnetic 
viscosity  or  diffusivity  compared  to  the  magnetic  viscosity  of  the 
quiescent  matter  above  and  below  the  zone.  Thus,  the  field  is 
confined  to  the  convective  zone,  resulting  in  traveling  dynamo 
waves  rather  than  a stationary  field  such  as  one  finds  in  the 
Earth  where  the  field  diffuses  out  of  the  core.  Prom  qualitative 
considerations  such  as  .illustrated  in  figure  2 or  from  the 
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analysis  in  sections  3 and  4*  we  find  that  the  waves  will  travel 

-I 

toward  the  equator  in  the  outer  hall  and  toward  the  pbles  in  the 
inner  half  of  the  solar  convective  zone.  The  velocity  of  propa- 
gation must  be  of  the  order  of  2 tn/ sec  in  order  to  give  the  ob- 
served 22 -year  magnetic  period. 

A wave  in  the  inner  half  will  decay  upon  reaching  the  pole 
due  to  the  lack  of  nonunifonn  I'otuLion.  it  will  introduce  only  a 
small  perturbation  into  the  outer  half,  which  will  be  amplified 
as  it  propagates  back  to  the  equator.  As  the  wave  nears  the  equa- 
tor it  will  decay  because  of  the  decrease  of  the  cyclonic  and 
anticyclonic  motions  associated  with  the  radial  convective  motions. 
The  wave  again  serves  only  to  introduce  a small  perturbation  down 
into  the  lower  half  of  the  zone,  which  is  then  amplified  as  it 
propagates  toward  the  pole. 

Initiating  the  wave  at  the  polar  end  of  the  outer  half 
or  the  equatorial  end  of  the  inner  half  of  the  convective  zone  is 
obviously  an  inefficient  process.  Hence,  considerable  amplifi- 
cation is  required  en  route.  Thus  we  conclude  from  (18)  that 

>>  (20) 

Neglecting  ^ k * in  ( 18 ) the  amplification  is  sx^s  [t  )k,  H V/2  i "]^ 
and  the  velocity  of  propagation  is  j HP  / (2  kj|  . Thus,  if 

t is  the  time  required  to  propagate  the  wavelength  2-n/k,  J we 


k,  HP 
2 


have 


1 h r i = 

1 2Z  1 

J 21  k,  1 

1 k.  1 

•x  In  figure  3,  for  instance , "T(  ■")  O , B„>  Gj  O in 

the  upper  half  of  the  convective  zone.  'Thus  HP  < O' and  the  wave 
propagates  in  the  positive  oc.  or  southerly  direction. 
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or 


t 


<i  H P I 
2 / 


Vk 


(21) 


The  amplification  is  o.p(Zr:)  . The  sunspots  and  other  solar  mag- 

P 

netic  activity  indicate  that  the  10  km  distance  from  the  polar  to 
the  equatorial  regions  is  about  one  wavelength,  27j/k,  . Thus, 

the  amplification  is  of  the  order  of  sxpfcny  or  500  in  each  half 
of  the  convective  zone.  Assuming  that  this  amplification  is 
sufficient  to  make  up  for  the  losses,  we  see  that  the  initial 
fields,  however  small,  will  be  built  up  after  a number  of  cycles 
to  where  they  limit  the  motions  resulting  from  the  Coriolis  forces 
and  reduce  jk,  HP  / 2 I relative  to  ^ * k A . When  this  point  is 
reached,  the  dynamo  henceforth  operates  in  a dynamically  stable 
state. 


If  more  were  known  concerning  the  dynamical  state  of  the 
convective  zone,  it  would  be  fruitful  to  use  turbulence  theory  to 
relate  HT  and  to  the  convective  velocities,  etc.  Suffice 

5 

it  to  say  that  if  the  depth  of  the  convective  zone  is  10  km, 
convective  velocities  of  the  order  of  1 m/sec  seem  to  be  ample. 
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6.  The  Dynamics  of  Hydromagnetic  Dynamos. 


Let  us  investigate  in  a semi-quantitative  manner  the 
dynamics  of  hydromagnetic  dynamos,  both  stationary  and  migratory. 
Consider  a fluid  body  of  magnetic  viscosity  v and  angular  velocity 
0.  Let  the  fluid  be  convecting  radially  with  a mean  velocity  v. 

We  denote  the  mean  toroidal  field  by  Bt  and  the  mean  poloidai 
field  by  . Let  the  scale  of  the  fields  be  L . Finally,  let  0 
be  the  angular  distance  from  the  north  pole. 

The  force  per  unit  volume  exerted  on  the  fluid  by  a mag- 
netic field  B is  (V  X B)  X • Thus,  the  force  exerted  by 

the  tordidal  field  against  the  nonuniform  rotation  drawing  out  the 
toroidal  field  from  the  poloidai  field  is  the  (p  component  of 
(Vx  Jg)  * *&/£**-  and  will  be  represented  by  (^.L)  . The 

force  exerted  by  the  magnetic  field  opposing  the  fonnation  rif 
a loop  from,  the  toroidal  field  will  be  represented  by  the  same 
expression,  Bt  3,.  LS)  * The  Coriolis  force  available  to  draw 
out  the  toroidal  field  is  v « si'N©  and  to  kink  the  toroidal 

field  to  form  loops  0)  where  a and  j3  are  dimension- 

less constants.  We  see  that  to  draw  out  the  toroidal  field  we 
must  have 

' ^ L 


to  produce  loops 


_ Bt  BP 

V /O  C o S 0 L 

Were  It  not  for  the  magnetic  field  inhibiting  the  Coriolis 
motions,  the  velocity  appearing  in  the  generating  term  Vx  (&  x t) 
would  be  proportional  to  the  Coriolis  forces  and  we  would  write 
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sin  Br 


OJ  V'/S  c OS 


®)  £* 


for  the  rate  of  generation  of  the  toroidal  and  poloidal  fields 
respectively!  b and  c are  constants.  However,  the  inhibition  of 
the  Coriolis  motions  by  the  magnetic  field  is  not  negligible, 
and  we  shall  therefore  consider  as  effective  only  that  fraction 
of  the  Coriolis  forces  which  is  in  excess  of  the  inhibiting  mag- 
netic forces.  The  generating  terms  of  the  toroidal  and  poloidal 
fields  become  _ 

i.  Bp.  [ u V 3 in  6 - J Bt  ^ 7^c  04  ® "■ 

The  dissipation  *-  V JB  may  be  writtin  as  e*et/L  and 
for  the  toroidal  and  poloidal  fields,  respectively. 

The  time -dependent  equations  for  Bt  and  are 


fct  L-1  ‘ 


L 

to 


f /f*-  L. 


> (22) 


~ 4-  ■ B>i  j v/3  cos© 

>v  < L1  r L-  L 


- 

L 


L 

Consider  first  the  equilibrium  values  of  3t  and  3P  , 
which  we  shall  denote  by  T and  TT  respectively.  Setting  the 
left-hand  sides  of  (22)  equal  to  zero  and  letting 

,,  = II  W = 1TT  . n‘-  ^ 


(23) 


r j-  ^ ^ w~) 

— -1-  | a,c  0's  S 4 a,  sir,  G rt  j c os  © — ^ J S'ia  0)  4-  4 ^ ^ ' J 

l^u  “\ 

— a,  Cos©  — si*  6 [(a,  COS©  - SieS)  4-  *4  J ^ 


(24) 


where 


, = ^<=>  o-i  V1 yiS  i_ ^ 3 ! — V cx  L-^-n, 

Let  us  assume  large  magnetic  Reynolds  numbers  so  that 


(25) 
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b.1,  a.'  >> 


ka  c 

Then  except  ir  the  neighborhood  of  the  special  value  of  & for 
which  a,  cos,  © = a*  ■air,  9 we  have 

(^a,co^8  — £vsir\©)  > ■>  u\ 

Expanding  (24)  in  powers  of  be  (a.co*©  - s 


we 


obtain 


w 


=.  a,  cas  !3  — 


l_l  =■ 


r1 


be.  (a,  «7r>  © - A ,e  OS©^  bvt“-  (»*s!n©  - S.tosOj 

•j  + O C « i j 


-a  +o  T^aO 


c^a^si'n©  - a,  cos©)  beh(ai  sin  © - a,  cos©) 

for  tslo©  > a,cosG  and 


J 


(26) 


vr 


= ^Sln© 


+■ 


v< 


bc.(a,cos<S  - 3j  &)  bvc  1 (3,e  os  © - av  *.>,©}* 


+0  ^3<,  av^ 


(27) 


u = 


_ _b_ 


(a. cos©  sin©)-h 


.V1 


—*^6*Ca,,aO 


c (ai  cess©  - Si  sin©)  be  ‘ (a  ,<-os©  -avsf  n wj* 
for  cos6  3i  sin  9.  The  corresponding  values  for  T and 

TT  may  be  competed; 

eV^Slvn©  > S i cot  0 • 

T*  =-£-  a.cMefa.sine-a.cosS)  P ^ /lA- 1 SfA^  ) + <M 

V b >v  ^,Sin8  - a,coa©y  v / 


. . *■ a i cos© 

c ( aisin © — • a,cos 


_ _ ^t*-**-'*  -a  i sim  6 « ✓ 

Vac.'  (mi  si<\  © - a ,c<as  9^a  'J 


; 28 ) 


rn  = yis,s'"8  ’ 


+ O(') 


J 


a , c OS  © > &*.  x(  r<  © ; 


-19- 


~TZ  9 5.,’a  G 


— ~ ros._€L. 


— +•  O (0,30 


bC“  ltos9“4»^s)  b'c  (a,  e os©  -a.**  ©/ 

H*- ^ a.s!r;6  (s,t6s0  •» 

c \ a,£o^-.aiSI„e/+U 


(29) 


x = f^±^y 


i»  + O"  " ( &,  j 3.) 


T~^  ' *“  ' (a,  coi6  '3\Si«S)> 

To  demonstrate  the  stability  of  these  equilibrium  states, 
consider  a small  perturbation  of  the  system  from  equilibrium.  We 
let  the  toroidal  field  be  represented  by  "T*  +• t and  the  poloidal 
field  by  TT 4-tc  . ~T  and  TT  sat  isfy  the  equilibrium  equations,  and 
z and  vc  are  assumed  to  be  small.  Substituting  into  the  time 
dependent  equations  (22)  we  obtain 


-&tt  _ 
•&»! 


= A.  K + A.vT:  +-  Ol(Tci‘t)  ; ^ = Ax.  * + AxxT  + O (30) 


where 


= - - .kl1  A.a  , (a.eos©  - 2TTl)^ 

’ LA  ^|_A  ' V " I 


A.,  = (»..-«  - 2Tn)  a..  = -fi  - -^4 

L.  " L /-**-  i— 


(31) 


As  a solution  we  let 


■n  - k,  OJ  t T = ka  «.XK3  UJ 


Then 


(32) 


(33) 


o = k,  (a,  - *-~0  + j o •—  l 1 A,,  f-  oj) 

which  yields  the  characteristic  equation 

oj  =-  j^(A.  + An)  i.  |\A..  - A + ^A.xAj.J  ^ >34) 

Computing  the  A < j from  (28),  (29),  and  (31)  we  obtain 
two  roots,  viz. 

(35) 


co-  = “ + O'  ( 3,3.) 
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cj>x  = 


be.  a,  cos  © ( Sk  V>n  0 — a,cos0y'  + 0(0 
uVL1 


for  ais'mS  ^ Si  cos  6)  and  for  a,  cos  © > a i *.  in  © 


ej 


= - + °_1  (s*,  s0 


(36) 


(37) 


St  Sil^S  (?S|CC5sS 


St  s ( h £*")  *f-  003 


(38) 


giving  an  exponential  return  to  the  equ 
(35)  and  (37)  represent  free  decay  back 
(38)  represent  active  degeneration  back 


iiuu  in  all. 
to  equilibrium; 
to  equilibrium. 


cases. 
(36) “and 
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7.  The  Dynamics  of  the  £>olar  Dynamo 


In  the  above  section  we  developed  in  a semi-quantitative 
manner  the  dynamics  of  hydromagnetic  dynamos.  The  calculations 
were  applicable  to  both  stationary  and  migratory  dynamos,  though 
our  main  interest  is  in  the  migratory  dynamo  at  the  moment.  From 
(28)  we  see  that  the  equilibrium  values  of  the  toroidal  and  pol- 
oidal  fields  may  be  written 


— p-2  _ !=■  — j V CQS0  +"  O ^ 

rr"  = J 

near  the  equator  where  a.'M*© >a^cK8  and  from  (29) 


(39) 


~T*‘  ^ O ‘ ^ 

■py  ' __  la  5>in  B COS  & +•  O ‘ J 


(40) 


near  the  poles  where  a,coa6>  > s.s.'nfi 

First,  we  note  that  we  do  not  have  the  purely  kinematical 
result  (Bullard.  1954)  that  there  is  a critical  velocity  of  the 
convective  motions,  which,  if  exceeded,  results  in  an  unlimited 
exponential  growth  of  the  field.  The  reason  is  that  the  dynamo 
results  from  Coriolis  velocities  and  not  from  the  primary  radial 
convective  velocities?  The  convective  velocities  give  rise  to 
Coriolis  forces,  resulting  in  a dynamo  only  insofar  as  the 
magnetic  field  does  not  prevent  the  manifestation  of  the  Coriolis 
forces  as  cyclonic  velocities  and  nonuniform  rotation. 

Second,  supposing  that  /Oc~>v  is  large,  we  see  that  near 
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the  poles  equilibrium  occurs  for  T~l  I . whereas  near  the 

equator  ~T~  '>s>  7T  0 The  crossover  occurs  where  s,co:>©-2jWbr , if 
a = (3,  at  9 = 45°.  ^hus,  above  middle  solar  latitudes  we 
expect  the  magnetic  activity  to  arise  primarily  from  poloidal 
fields  and  below  middle  latitudes  from  toroidal  fields# 

Let  us  organize  the  above  dynamical  considerations  into  a 
coherent  picture  of  the  solar  dynamo#  The  dvnamo  was  discussed 
in  section  5 from  a purely  kin9matical  point  of  view:  In  the  upper 

half  of  the  solar  convective  zone  we  have  a relatively  weak  travel- 
ing dynamo  wave  originating  near  the  poles  and  being  strongly  am- 
plified as  it  travels  toward  the  equator.  (39)  and  (40)  imply 
that  it  will  first  be  the  poloidal  field  that  is  the  more  strongly 
amplified*  somewhere  in  the  middle  latitudes  the  strong  amplifi- 
cation of  the  poloidal  field  will  cease  and  the  toroidal  field  will 
be  amplified.  We  expect  then  that  magnetic  phenomena  arising  from 
the  poloidal  f ield  will  be  observed  primarily  toward  the  pole 
from  middle  latitudes  (though  of  course  the  poloidal  field  may  be 
strong  enough  to  have  some  observable  effects  nearer  the  equator); 
phenomena  arising  from  the  toroidal  field  will  be  observed  toward 
the  equator  from  the  middle  latitudes#  Pew,  if  any  local  magnetic 
phenomena  such  as  sunspots,  flares,  and  prominences  should  be 
observed  at  the  equator,  where  both  the  poloidal  and  toroidal 
fields  will  be  decaying# 

Now  the  observed  nonuniform  rotation  of  the  sun,  wherein  the 
surface  equatorial  regions  rotate  more  rapidly  than  the  average  cf 
the  sun,  ( contrary  to  the  effect  produced  by  the  Coriolis  forces 
of  the  convective  motions),  can  be  interpreted  in  the  3ame  manner 


as  Bullard  et  al,(1950)  has  interpreted  the  corresponding 
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terrestrial  phenomenon,  the  westerly  drift  of  the  Earth’s  magnetic 
field.  The  nonuniform  rotation  depends  on  the  strength  of  the 
toroidal  field,  and  hence  theoretically  should  increase  toward  the 
equator  from  the  middle  latitudes  where  it  is  essentially  zero. 

This  is  as  observed  (Kuiper,  1953). 

The  northern  and  southern  hemispheres  of  the  sun  operate 
as  independent  dynamos  except  for  whatever  loose  couplir.g  ray 
exist  due  to  fields  diffusing  across  the  equator.  The  opposite 
polarity  of  sunspots  in  the  two  hemispheres  indicates  that  the 
coupling  is  such  that  the  toroidal  field  changes  in  sign  across  the 
equator.  The  looseness  of  the  coupling  between  hemispheres  and 
between  the  dynamos  of  the  Inner  and  outer  halves  of  the  convec- 
tive zone  explains  the  lack  of  complete  similarity  between  hemi- 
spheres and  between  consecutive  cycles, 

* 

Elsasser  has  suggested  that  sunspots  arise  as  a cotr 
sequence  of  a strand  of  the  overall  mat  of  toroidal  field  rising 
through  the  photosphere.  This  gives  an  inverted  ”UW  shaped  loop 
of  flux  rising  from  the  main  toroidal  field  up  to  the  surface  of 
the  sun,  and,  together  with  the  migratory  dynamo  yielding  magnetic 
waves  of  alternating  sign  traveling  from  the  poles  to  the 
equator,  gives  a sunspot  cycle  with  the  proper  reversal  of  spot 
polarity  between  the  half  cycles  of  11  years. 

The  question  now  arises  as  to  why  we  should  obtain  a 

% 

migratory  dynamo  for  the  sun  and  a stationary  dynamo  for  the 
core  of  the  Earth.  First  we  note  that  below  and  above  the  con- 
vective zone  there  is  presumably  not  sufficient  turbulence  to 
significantly  diffuse  magnetic  fields.  Beneath  the  convective 

^Unpublished 


-25 


zone  the  conductivity  is  of' the  order  of  10  ’’  mho/m  (Cowling,  1945 )». 

The  magnetic  viscosity  becomes  10”^  m^/sec  and  the  time 

K 9 

required  to  diffuse  oven  10°  km  is  of  the  order  of  10  years* 

The  conductivity  above  the  convective  zone  is  also  large,  and  dif- 
fusion upward  is  ruled,  out  unless,  due  to  some  local  anomaly,  the 
field  can  poke  out  through  the  less  dense  overlying  layers  pushing 
tha  matter  aside,  as  evidently  happens  locally  in  sunspots  and 
prominences.  We  conclude  that  the  magnetic  fields  generated  by 
convection  will  be  confined  largely  within  the  convective  shell 
itself. 

With  the  great  length  of  the  convective  zone  as  compared 
to  its  depth,  we  do  not  expect  the  magnetic  activity,  whose  scale 
is  determined  by  the  relatively  shallow  depth  of  the  zone,  to  be 
correlated  over  the  entire  solar  sphere.  Hence,  a magnetically 
coherent  region  of  the  sun,  e.g.  the  northern  solar  temperate 
zone,  is  isolated  and  is  in  an  effectively  unclosed  space.  This 
is  in  contrast  to  the  situation  in  the  core  of  the  Earth,  where  the 
depth  and  circumference  of  the  core  are  of  the  same  order  of 
magnitude,  resulting  in  magnetic  coherence  over  the  entire  closed 
space.  It  is  the  solution  of  (4)  and  (7)  in  a closed  or  unclosed 
space  that  yields  a stationary  or  migratory  dynamo  respectively. 
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8.  Magnetic  Stars 

Babcock  (1953)  has  shown  that  many  A and  M stars  have 
general  magnetic  fields  of  several  thousand  gauss.  Several  cases 
are  known  to  be  magnetic  variables  with  periods  from  several  days 
up  to  a year.  Babcock  states  that  in  no  case  is  there  good  evidence 
that  the  magnetic  field  is  without  time  fluctuations,  though  many 
stars  exhibit  overall  magnetic  fluctuations  which  are  not  periodic. 
Given  a star  with  a deeper  and  more  active  convective  zone  than 
the  sun  there  seems  no  basic  difficulty  in  accounting  for  the 

i 

extreme  magnetic  activity  that  Babcock  observes.  Prom  (39)  and 

(40)  we  see  that  the  field  is  proportional  to  L/=>  ^ v • If 

the  eddy  diffusivity  v is  proportional  to  Lv  we  have  L ^wv'/l, 

Increasing  L. , and  v by  only  a factor  of  ten  from  the  10  km. 

5 x -10”  sec”  , and  1 m/see  used  in  the  sun  gives  fields  with 
2 

intensities  10  time3  the  solar  fields  and  magnetic  periods  of  two 
yearsj  increasing  v to  1 km/sec  gives  fields  10  times  the  solar 
field  with  a period  of  one  week. 

Many  magnetic  stars  exhibit  a marked  asymmetry  in  their 
variations.  Indicating  that  they  have  a large  stationary  component 
of  their  magnetic  fields.  The  discussion  at  the  end  of  the  previous 
section  would  imply  that  their  convective  zones  are  proportionately 
deeper  than  in  the  sun,  resulting  in  a large  stationary  dynamo 
component  along  with  the  migratory  dynamo. 
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Appendix 


To  consider  further  solutions  of  (13)  we  let 

t5  (fej  t ) = p Cfcj 

Substituting  into  (13)  we  obtain  for  the  equation  for  P (fc> 


(41) 


( L - 


(42) 


« ?-(k'-k')  p ik'; 

® V-  Ck 

Thd  coupling  between  the  modes  and  J5*  is  given  by 

ii/,  _ (-  e»-  k*  + g»  0 * JL  (k  - ^ ')  ^ Ck  - A*}  (43) 

' (wi  *J  tlX  ^ f ^ k'  <) 

u (fc,  fc7  is  a monotonically  and  rapidly  decreasing  function  of 

both  k and  k'  because  of  the  two  factors  in  the  denominator. 

Usually  |h  and  are  rapidly  decreasing 

functions  of  the  magnitude  of  their  arguments.  Thus  U (§»_,  K>  Will 

be  small -unless  |P — jk  ” . This  means  that  the  coupling  to  a 

given  mode  Jk  usually  does  not  extend  over  a large  neighborhood 

of  jk  „ which  allows  an  approximate  solution  of  (42)  without  undue 

c alculation. 


(43) 


Consider  the  solution  of  (42)  when 

v*  = s,  h , ret)  - r e:  (*; 

where  H and  P are  constants.  Their  Fourier  transforms  are 

k (t;  = H s (t; , j (k)  = r & oo  e.  oo 

where 


(44) 


(45) 
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(42)  reduces  to  the  linear  homogeneous  integral  equation 

r 

p — X(l'-i)  J dkx  PC  ( kL  J p (k0  k»/)  (47) 


x uj  = £ h r 


/ i i / ) 

I k — ' L l / 'N  _ . ^ ^ k.  ^ ~ k»  .z 


(48) 


(47)  has  nontrivial  solutions  (Lovitt,  1950)  if 

v /,  X rlT  k"^L  , \ j X*(k,)  lK(^,lzJK(ki;ki)  | 

°-  1 - *-  ~rrl  j^/K  (u;.QKa-t>oi 

-r  * * • (49) 

(48)  is  an  absolutely  and  permanently  converging  series  in  Xlk^ 
Thus  to  a first  approximation  we  drop  all  terms  second  order  and 
higher  in  X (.k,)  and  obtain  the  eigenvalue 

x(t.)  = [ ^f«U*  K .ujj)  (50) 

Using  (48)  and  (50)  we  obtain 

/ / i ».  \ 3^*- 

Hr  - - i + ( 51) 

x e(o)  k, 

Solving  for  (0  we  obtain 


(51) 


tc  e(Q)  k, 


- [i  ^22 
1 ' L 5v 


3) k,  Hr 
vV 


(52) 


With  X (.k, ) given  by  (50)  the  solution  of  (47)  is  given 
by  the  absolutely  and  uniformly  converging  s eries 

/"  I K k.)  K(k>,  k) 

P ^ w A(OK(^k.)-  ko  Mk|k)/  '5'5 

for  suitable  k.  . We  shall  put  k«  = O . Then,  to  a first 
approximation 

P (fc)  ~ XCkXCk,0)  = i HP  j ^ Hey 


(54) 
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Now  iVl  - - I 


or  j?  V3  , The  first  case  gives 


a degenerating  traveling  dynamo  wave  and  the  solution  vanishes  for 
large  t • Putting  i3/i  — — we  obtain  the  regenerating 


traveling  dynamo  wave 

b(sk,o  = i ht 


ki  £, 


(|  * 


L & J i 1 


* cxP^/--( ---k'  — ) U e*f= 


£ V** 


J 


(55) 


traveling  with  a speed 

V3  / Tt  g co;  k , H r ^ 

\ 2 y 

in  the  direction. 

In  a similar  manner  we  may  solve  (31)  in  the  cas; 

Vv  = 6,h  jco  , rcE;  = r . 

Their  Fourier  transforms  are 

b Cjfc)  = H SC fc.)jO0  , ^ (V  = P ^(feJ 


( 56) 


that 


ivVor»ft 


j (k.)  = -5-  d H J (2)  . 

«»  — Ota 

(42)  again  reduces  to  a linear  homogeneous  integral 
equation.  Defining 


L ) — 


C^+  ^k^x^^k'x) 


we  obtain 


= X y=^^i  L.  Lj  Lv  ) p>  ( k,j  ) . 


To  a first  aonroximation 


X ( k')  — J dkj  L (,1 k-Oj 


- /•  ■ 

_ c ( ^ ^ fc»  / 

?r  j co) 


and 


*/* 


'k  ‘)co^  h r 

3 -yV 


To  a first  approximation 

, ^ . | j r->  l< . j 0 

p ^ - ■ H 1 ~T^r^) 

Neglecting  the  degenerative  solution  resulting  from  i */s  = - I J 
we  have  the  regenerative  traveling  dynamo  wave 


X «t3t_p 


■y  ) < o)  k , H 
2 V*1 
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